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The capacity expansion planning problem of the renewable energy industry involves decisions 
regarding the optimal mix of different plant types, locations where each plant should be built, and 
capacity expansion decisions over the planning horizon for each plant. The aim of this paper is to 
develop a goal programming model, based on a multi-source multi-sink network, in order to locate five 
renewable energy plants for electric generation in five places located in the autonomous region of 
Cantabria, in the north of Spain. As different types of plants can be placed in each location, the goal is to 
locate one plant in each place, maximizing the number of plants that are matched with comparable 
locations, in a way that the total deviations from goals are minimized. 
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1. Introduction 


Mitigation of climate change, in order to reach the targets set 
by the Kyoto Protocol, requires significant reductions of Green- 
house Gas emissions and, consequently, the adequate policy 
formulation for fossil fuel energy substitution by renewable 
energy sources. In modern technologically developed societies, 
where the complexity of the interactions between the economic 
activities and environmental systems are recognized, strategic 
decisions must be made under an increasingly complex and 
unstable environment, characterized by a fast pace of technolo- 
gical evolution, changes in market structures, and new societal 
concerns [1]. These problems demand the use of analytical tools 
that describe and evaluate the problem in its social, environ- 
mental, economic and technological dimension. 

The capacity expansion planning problem of the renewable 
energy industry is not only concerned with the selection of the best 
renewable energy project among different alternatives, but also with 
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developing expansion plans for the generation, transmission and 
distribution systems. Such plans involve decisions regarding the 
optimal mix of different plant types, locations where each plant 
should be built, and capacity expansion decisions over the planning 
horizon for each plant. A multitude of technical, financial, environ- 
mental, legal, social and political objectives and/or constraints, some 
of which are not even quantifiable, some of which are conflicting, 
participate in the decision process, rendering capacity expansion 
planning in the power field a very complex issue [2]. 

In this paper a Goal Programming model, based on a multi- 
source multi-sink network, is presented. The paper is organized as 
follows. In the next section the Goal Programming methodology is 
shown. Next, the model is applied to locate five renewable energy 
plants for electric generation in five places located in the autono- 
mous region of Cantabria, in the north of Spain. Finally, a concluding 
section with the main results of the paper is presented. 


2. The goal programming model 


Although Goal programming is itself a development of the 
1950s, it has only since the mid-1970s that Goal programming 
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received truly substantial attention. Much of the reason for such 
interest results from Goal programming’s demonstrated ability to 
serve as an efficient and effective tool for the modeling, solution, 
and analysis of mathematical models that involve conflicting 
goals and objectives [3]. 

Goal programming can be used to help satisfying decision- 
makers to attack problems in which a multitude of conflicting goals 
exist. During recent years, the increasingly complex social, eco- 
nomic, technological and environmental scenario has made multi- 
criteria analysis a valuable tool in the selection of renewable energy 
projects. As a result, multi-criteria decision-making methods such as 
Analytical Hierarchy Process, ELECTRE, PROMETHEE, TOPSIS, VIKOR, 
have been extensively used in the selection of renewable energy 
projects in areas such as wind-farm projects, geothermal projects, 
hydro-site selection, etc. [4-14]. 

Whereas these methods have been widely used in the renew- 
able energy industry, the use of the Goal Programming metho- 
dology in the selection of renewable energy plants, locations, or 
capacity expansion decisions has been limited to a few works. 
Koroneos et al. [15] applied a multi-objective model on the island 
of Lesvos, Greece, where various renewable energy sources can be 
exploited with the existing conventional systems to satisfy the 
needs of the island’s economy. Daim et al. [16] developed a fuzzy 
goal programming model for the state of Oregon to create a 
renewable energy portfolio with the objective of responding to a 
25 percent of the electricity demand by renewable resources in 
2025. Jinturkar and Deshmukh [17] developed a fuzzy mixed 
integer goal programming model for rural cooking and heating 
energy planning in Central India. 

In linear programming, only one objective is permitted, thus, 
even though multiple goals may confront the decision-maker, all 
progress towards these goals must be measured on a common 
scale, often profit or cost. The largest drawback of linear pro- 
gramming is that it requires unidimensionality in the objective 
function; that is, all goals must be expressed in common units and 
combined to give an overall single measure of effectiveness. The 
requirements are that the decision-maker must be able to estab- 
lish and prioritize the goals and to express the relationships 
between the decision variables and goals with linear functions. 
Once the decision-maker has provided an ordinal ranking of the 
goals, the goal programming model minimizes the deviations 
from the goals subject to the constraints that have been set. 

The general form of a Goal Programming problem may be 
expressed as 


Minimize X (df +d; ) 
i 


variables in the goal equations, b(i=1,2,...,m) are the targets or 
goals, and a, are the coefficients of the variables. 


3. Application 


The aim of this section is to develop a goal programming 
model in order to locate five renewable energy plants for electric 
generation in five places located in the autonomous region of 
Cantabria, in the north of Spain. Seven plants (alternatives) are 
candidate to be placed in these seven locations: La Braguia, La 
Vega, Estacas, Lunada, Potes, Reinosa, and Santander. The attri- 
butes considered in order to evaluate the designed systems are: 
Power (Gw), Investment (J), Tons of emissions of COz avoided per 
year (TCO3/y), Operation and Maintenance costs (OM), Jobs (J), 
Distance between plants (D), and Social acceptance (S). Table 1 
shows the data for the first five criteria [18], Table 2 shows the 
distance between locations and Table 3 the social acceptance 
expressing the overview of opinions related to the energy systems 
by the local population regarding the hypothesized realization of 
the renewable energy projects. Let us consider for analytical 
purposes, that the index social acceptability has been expressed 
by the local population in a scale of 1 (low acceptance) to 10 (high 
acceptance). 

Table 4 shows the comparability between plants and locations. 
As different types of plants can be placed in each location, the goal 
is to locate one plant in each place in a way that the total 
deviations from goals are minimized. Let us consider that the 
autonomous government establishes the following goals: the 
power generated, emissions avoided and jobs created must be 
higher than 110 x 10° Gw, 18 x 10° t, and 70 jobs; with regard to 
investment and operation and maintenance costs, these figures 
must be limited to 100 x 10° € and 350,000 €/year respectively. 
The distance between plants must be maximized to 2910 Km and 
the social acceptance must be as close as possible to the highest 
level (50). 

Fig. 1 shows the appropriate network, where each plant is a 
source node and each location is a sink node. By adding a super 
source node connected to each source and a supper sink node 
connected to each sink, the multi-source multi-sink problem is 
transformed into a maximum flow problem that maximizes the 
number of plants that are matched with comparable locations. 


Table 2 
Distances between locations (Km). 


La Braguia La Vega Estacas Lunada Potes Reinosa Santander 


s.t. 
n P _ La Braguia - 6.5 15 39 126 61 53 
X (ajx)d; +d; =b; La Vega = 8 32 121 54 65 
j=l Estacas - 24 128 62 72 
ss Lunada - 143 78 56 
+ 
Xj+d; d; =0 0) Potes - 130 107 
where d; is the amount by which goal i is underachieved, d} is cana 7 13 
the amount by which goal i is overachieved, x;(j=1,2,...,n) are the 
Table 1 
Alternatives for electric generation [18]. 
Alternative Gw I TCO3/y J OM 
Wind power (10 < P < 50 MW) 60,959,000 23,425,000 9,649,680 15 37,750 
Hydroelectric (P < 10 MW) 3,940,100 7,500,000 472,812 8 7250 
Hydroelectric (10 < P < 25 MW) 962,000 14,000,000 255,490 8 14,000 
Hydroelectric (25 < P < 50P MW) 412,000 21,035,000 255,490 12 21,000 
Solar (P> 10 MW) 306,031 5,320,000 482,856 10 1792 
Biomass (P <5 MW) 13,612,500 9,015,000 2,524,643 15 27,100 
Biomass (P > 50 MW) 37,770,000 47,958,400 4,839,548 20 75,000 
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There is an arc joining the super source to each plant, and an arc Fso,1—F1,3—F1,10—Fi,11 = 0 (17) 
joining each pair of compatible matches, and an arc joining each 

place to the super sink. Because arcs do not exist between Fso,2—F2,8—F2,9—F2,12 = 0 (18) 
noncomparable mates, a flow of k units from the super source 

to the super sink represents an assignment of plants to places in Alternatives Places 

which k compatible couples are created. Because the arc joining 

each place to the super sink has a capacity of 1, conservation of o 

flow ensures that each place will be matched with at most one 

plant. Similarly, because each arc from the source to a plant has a \ 

capacity 1, each plant can be pared with at most one place. The (2) ASNS Cc 

formulation for this goal programming model can be expressed as Se 


Minimize J (di +d} +d; +d, +d} +d5 +d7) (2) SOK? 
i VAX X 


KD 
s.t. 7 
60,959,000F1 8+... +37,770,000F7 14+dī—d] =110%10° (3) (+) RIT 


\y 
KS Sink 
23,425,000F),¢+...+47,958,000F7,14+d5—d} = 100%10° (4) K AK 


XAN 
9,649,680F1 8 +. . .+4,839,548F7 14+d3 -df = 18%10° (5) TAKE ee 


\ 
15Fi.g+...+20F714+d,—-d}{ =70 (6) Ys 
18+... 714+d,—-d4 = 
O222°6 


37,750F1 8+.. .+75,000F7 14+dz -dł = 350,000 (7) 
6,5Fgsit...+74Fy4si+ds—d¢ = 2910 (8) O) © 
7Fı s+... +9F714+d; -d7 =50 (9) 
Fso si—F. F F F. F F Fso7 =0 10 
BOSE SO; i i E 50,6490 (10) Fig. 1. Multi-source multi-sink network. 
Fg oj +Fo,si + F10,si +F11,si + F12,si + Fsi,so = 0 (11) 
F1,8 +F2,8 +F3,3 + Fa,s—Fs,si = 0 (ia. Tables 
Places to locate the plants. 
F29 +F4,9 +F6,9—F9,si = 0 a3) Alternative Location 
F110 +F3,10 +F5,10 +F6,10—F10,si = 0 a4 (1) Wind power (10 < P < 50 MW) Estacas (10) 
(3) Hydroelectric (10 < P <25 MW) Santander (14) 
F111 +F3,11 +F5,11 +F6,11 +F7,11—F11,si = 0 (15) (5) Solar (P> 10 MW) Potes (12) 
(6) Biomass (P < 5 MW) Reinosa (13) 
7) Biomass (P > 50 MW Lunada (11 
F212 +F5,12 +F712—F12,si = 0 (16) (7) ( ) un 
Table 3 
Social acceptance (1-10) of each plant in each location. 
Alternative La Braguia (8) La Vega (9) Estacas (10) Lunada (11) Potes (12) Reinosa (13) Santander (14) 
(1) Wind power (10 < P < 50 MW) 7 9 7 9 
(2) Hydroelectric (P < 10 MW) 5 6 5 7 
(3) Hydroelectric (10 < P< 25 MW) 6 7 6 7 
(4) Hydroelectric (25 < P< 50P MW) 4 5 6 
(5) Solar (P> 10 MW) 8 7 9 
(6) Biomass (P< 5 MW) 9 8 8 9 
(7) Biomass (P > 50 MW) 9 8 9 
Table 4 
Comparability between plants and locations. 
Alternative La Braguia (8) La Vega (9) Estacas (10) Lunada (11) Potes (12) Lunada (13) Potes (14) 
(1) Wind power (10 < P < 50 MW) * * * * 
(2) Hydroelectric (P < 10 MW) * * + * 
(3) Hydroelectric (10 < P < 25 MW) * * * * 
(4) Hydroelectric (25 <P<50P MW) * * * 
(5) Solar (P> 10 MW) * * 
(6) Biomass (P< 5 MW) * * * * 
(7) Biomass (P > 50 MW) * * + 
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Table 6 
Deviation from goals. 


Gw I TCO3/y J OM Distance Social acceptance 
Goals 110 x 10° 100 x 10° 18 x 10° 70 350,000 2910 50 
Results 113,609,530 99,718,400 17,752,217 68 155,642 2323 43 
d? 3,609,530 0 0 0 0 0 0 
dy 0 281,600 247,783 2 194,358 587 ri 
Fso,3—F3,8—F3,10—F3,11 = 0 (19) of plants that are matched with comparable locations, in a way 
that the total deviations from goals are minimized. The model 
Fso,4—F4,8—F4,9 = 0 (20) shows as Goal programming can be used as an efficient and 
effective tool to help decision-makers to attack problems in 
Fso,5—Fs,10—Fs,11—Fs,12 = 0 (21) d f . . 
eveloping expansion plans for the renewable energy industry, 
Fso,6—F6,9—F6,10—F6,11 = 0 22) where a multitude of conflicting goals exist. 
Fso,7—F7,11 -F7,12 = 0 (23) 
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where Eqs. (3)-(9) correspond to the power, investment, emis- 
sions, job, operation and maintenance, distance and social accep- 
tance goals; Eqs. (10)-(23) are the conservation of flow constraints 
for each node; Eq. (24) ensures that each place will be matched 
with one plant; Eq. (25) ensures that five plants must be located 
five locations; and Eq. (26) indicates that the flow variables are 
binary. 

The problem was solved using Lingo and the solution is 
given in Tables 5 and 6. Table 5 shows the places where each 
plant should be located and Table 6 shows the deviations from 
the goals. Goals 1, 2, and 5 are achieved, while goals 3, 4, 6 and 
7 are not achieved. The power generated is higher than the 
established goal, whereas the investment and operation and 
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4. Conclusions 


Mitigation of climate change requires the adequate policy 
formulation for fossil fuel energy substitution by renewable 
energy sources. In modern technologically developed societies, 
where the complexity of the interactions between the economic 
activities and and environmental systems are recognized, the 
capacity expansion planning problem of the renewable energy 
industry involve decisions regarding the optimal mix of different 
plant types, locations where each plant should be built, and 
capacity expansion decisions over the planning horizon for 
each plant. 

In this paper a Goal Programming model, based on a multi- 
source multi-sink network, is applied to locate five renewable 
energy plants for electric generation in five places located in 
the autonomous region of Cantabria, in the north of Spain. 
As different types of plants can be placed in each location, the 
goal is to locate one plant in each place, maximizing the number 
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